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Abstract 

Patterns of shape-phase transition in the proton-neutron coupled systems are studied within the 
5D-pair shell model. The results show that some transitional patterns in the S'-D-pair shell model 
are similar to the i7(5) — ^[/(S), f/(5) — SO{Q) transitions with signatures of the critical point 
symmetry of the interacting boson model. 
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I. INTRODUCTION 



Recently, based on the Generalized Wick Theory a nucleon-pair shell model(NPSM) 
was proposed[2|, in which nucleon pairs with various angular momenta are used as building 
blocks. Since modern computers fail for the calculation in the full shell model space for the 
medium weight and heavy nuclei, some truncation scheme need to be used. 

The tremendous success of the interacting boson model (IBM) 3|, suggests that S and D 



pairs play a dominant role in the spectroscopy of low-lying nuclear modes . Therefore, 

one normally truncates the full shell-model space to the collective S-D pair subspace in the 
NPSM. The latter is called the SD-paii shell modeKSDPSM)^, H, fl. 

A crucial point in the SDPSM is the validity of the S-D pair truncation. In Ref.j^, shell 
model foundations of the IBM was reviewed by lachello, the results seem to indicate that 
the S-D pair truncation is a reasonable approximation to the full shell model space. This 



problem was also studied in 



11 



12| with the conclusion that the S-D pair subspace works 



well in the vibrationa 
But Dr. Zhao's work 



re gion , but in the deformed region, the inclusion of G pairs is crucial. 



13 



14j show that the essential properties within the full shell model 



space survive in the S-D pair subspace. What's more, if a pure quadrupole-quadrupole 
interaction and a reasonable collective S-D-pair were used, the rotational behavior can be 
produced very well within the S-D pair subspace. The fact that the SDPSM can describe 
the collectivity of low- lying states for nuclei around A=130 15, 16, 17, 1^, 19, 20, 21 1 also 



imply that the S-D pair truncation is a good approximation to the full shell model space. 

Nuclei, as a mesoscopic system, have been found to possess interesting geometric shapes, 
such as spherical (vibrational (f/(5)), axially deformed {SU{3)), and 7-soft (0(6)), which 



is usually described in terms of the Casten triangle [22|. The search for signatures of tran- 
sitions among various shapes (phases) of atomic nuclei is an interesting subject in nuclear 
structure theory. An understanding of such shape (phase) transitions may provide insight 
into quantum phase transitions in other mesoscopic systems j^. 

Theoretical study of shape phase transitions and critical point symmetries in nuclei was 



mainly carried out 
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40| in 



the interacting boson model for identical system(IBM-I) 3|]. The investigations on nuclear 



shape phase transition and critical point symmetry for identica 



been carried out with fermionic degrees of freedom in 41 



42 



nucleo n sy stem have also 



43 



44 



45 



46|. Recently, 
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investigations o 
also carried out 



le shape phase transitions and critical point symmetries in nuclei were 



47 



481 ] in the proton-neutron interaction boson model(IBM-II) j 
Since the SDPSM is also built up from SD pairs, it is expected that the SDPSM can 
produce similar results to those of the IBM. Our previous work show that the vibrational, 



rotational, and gamma-soft spectra can be well reproduced (49| similar to the f/(5), SU{3) 
and 5*0(6) limiting spectra in the IBM. What's more, the vibrational-rotational phase tran- 
sition for identical system can also be produced within the framework of the SDPSM with 
fermionic degrees of freedom [s^. Since nuclei are neutron-proton coupled systems, and a 
rich phase structure can be obtained in the neutron-proton IBM, it is interesting to see if 
the phase transitional patterns in the neutron-proton coupled system can be produced in 
the SDPSM with fermionic degrees of freedom. This is the aim of this paper. 

II. MODEL 

In the shell model description, the pairing and quadrupole-quadrupole interactions are the 
most important short-range and long-range correlations. Considering that the Hamiltonian 
used to study the shape phase transition in the IBM is mainly composed of the monopole 



pairing and quadrupole-quadrupole interaction(e.g., Ref . |47l. l48l| ) . a schematic Hamiltonian 
is adopted in the SDPSM, which is a combination of the monopole pairing interaction and 
quadrupole-quadrupole interaction with 

a 



^11 'rij 



where X in Hx is denoted as f/(5), SU{3), or 5*0(6) corresponding vibrational, ro- 
tational, or gamma-soft limiting case in the model, G„ and k„ are the pairing and 
quadrupole-quadrupole interaction strength between identical-nucleons, respectively, k is 
the quadrupole-quadrupole interaction strength between proton and neutrons. In this pa- 
per, we set Gj, = and = k^. 

To study the phase transitional patterns, the Hamiltonian for proton-neutron coupled 
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system is written as 

H={l-a)Hui^^) + aHx. (2) 

where < a < 1 is a control parameter, Hx is taken as Hsu{3) when we study vibration- 
rotation transitional patterns, and is taken as Hso{6) when we study vibration to 7-soft 
transitional patterns. 

The E2 transition operator adopted is 

T{E2) = e.Q(2)+e.Q(2), (3) 

where e.j^{eiy) is the effective charge for proton(neutron). 
The collective S'-pair is defined as 

5t = ^y(aaO)(CtxCt)0 (4) 

a 

In this paper, the S'-pair structure coefficient, as an approximation, is fixed to be y{aaO) = 
^\/ n ^jv "' 'where fia is defined as fia = a + 1/2 and is the number of pairs for like-nucleons. 
The D-pair is obtained by using commutator 

= i[Q(2),^t] = J2y(ab2) (cl X Ciy. (5) 

ab 

After symmetrization, it is easy to obtain that 



y{ab2) = -ig(a62) 



The details of the model can be found in 



III. RESULTS 



y{aaO) y{bbO) 



(6) 



To identify shape phase transitions and determine the corresponding patterns, lachello et 
al initialed a study on effective order parameters, which should display different critical be- 
haviors for the phase transitions with different order. Specifically, the quantities related with 
isomer shifts, defined as V2 = (< Oj^ln^lO^ > — < Oj^ln^lO]^ >)/^ and vi, = (< 2f\h(i\2f > 



< 0^|firf|0^ >)/^ ) were proposed as effective-order parameters in 33|]. Consequently, 



some other quantities, such as the B{E2) ratios Ki = B{E2]A'[ 2|)/fi(^2;2| 0^) 



and K2 = B{E2; 0^ 2+)/5(E2; 2+ Of) |4oj as well as the energy ratio Rqo = E^+fE^ 
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were also suggested as the effective order parameters to identify phase transitions and the 
corresponding orders. Therefore, to study the shape phase transition in the S'-D-pair fermion 
model space, V2, v'21 in which the d-boson number operator is replaced by Z^-pair number 
operator N£, in the SDPSM, Ki, K2 and -Rgo be studied in this paper. Because the im- 
portance of -R42 = / Er^+ in determining the limiting cases and shape phase transitions 5l|] , 
i?42 is also presented. 



A. vibration-rotation transitional patterns 

We begin by considering the vibration-rotation phase transition. A system with = 
Ny = 3 in gds shell was studied. By fitting i?42 = E^+/E2+ = 2 for vibrational case, the 
parameters used to produce the vibrational spectra were obtained, and pre sented in Tabled 
The detailed discussion about the vibrational spectra can be found in |49f|. In the SDPSM, 
the full shell model space was truncated to the SD-paii subspace. The investigation on the 
validity of the S-D pair truncation in 1^, [ll, 12 1 show that the S-D pair truncation can 

nn n 

not produce the rotational spectra. But Dr. Zhao's work|13l. |14| | and our previous work[49| 
show that if a pure quandrupole-quadrupole interaction and a reasonable collective S-D 
pair were used, the rotational behaviors can be produced very well. It is found that with 
2k^ = 2k^ = k = 0.2MeV/r^, the similar resuhs as the SU{3)^ x SU{3)^ limit of the IBM 
can be produced, the typical energy ratios E^+ / E2+ and Eq+ / E2+ are 3.33 and 6.96, close 
to the IBM result 3.33 and 7. The detailed discussion can be found in Ref. [4^. 

TABLE I: The parameters used to produce the vibrational, rotational and 7-soft spectra. is in 
unit of MeV, Ko- and k are in unit of MeV/rg. 





limit 










K 


vibration-rotation 


vibration 


0.5 


0.5 








0.01 




rotation 








0.1 


0.1 


0.2 


vibration-7-soft 


vibration 


0.5 


0.5 








-0.01 




7-soft 


0.15 


0.15 








-0.015 



Energy ratios R42 and Rqq against control parameter a are shown in FigUl Fig{T^ shows 
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that the energy ratio R42 is 2(when a = 0) and 3.3 (when a = 1), typical values of vibrational 
and rotational spectra in the IBM 3(. It is also shown that the rapid change occurs when 
0.3 < a < 0.6, which indicates a phase transition occurs in this region. 

The energy ratio Rqq given in Fig{T)3 shows that similar behavior as that of the IBM for 
finite number of boson Nb is reproduced. It exhibits a modest peak followed by a sharp 
decrease across the phase transition, a typical signature of the Ist-order quantum phase 



transition 



52|. 



The SDPSM results of V2, V2, Ki and K2 are given in Figl2] and Fig|3l The effective 
charges were fixed with = 3ejy = 1.5e. As argued in |33|, U2, v'2 should have wiggling 
behaviors in the region of the critical point due to the switching of the two coexisting phases 
for the first order phase transition, then the obvious wiggling behaviors shown by V21 in 
Figl2]further confirm the transition is first order. The results of B{E2) ratio Ki is consistent 
with those of other effective quantities {33!. I40I . The critical behavior of K2 seems to deviate 
from the character of the first order phase transition. 

In the IBM, the critical point symmetry {29! between f/(5) and SU{3) is X{5). Since 
the shape phase transition between vibrational and rotational limit can be reproduced in 
the SDPSM, it is interesting to see if the properties of the X(5)-like symmetry also occurs 
within the SDPSM. We found that there is indeed a signature with a = 0.54 in the SDPSM 
similar to that of the X{5) in the IBM. A few typical values are given in Tab le ITTl from 
which one can see that typical feature of the X{5) symmetry stated in Ref. 52, l53|] indeed 
occurs in the SDPSM. For example, R42, Rm and Eq+ / E2+ is 2.91, 1.05 and 5.32 in the 
SDPSM calculation, close to the IBM results 2.91, 1.0 and 5.67, respectively. 



B. vibration-7-soft transitional patterns 

The investigation on vibration-7-soft shape phase transition in the IBM has been studied 



m 



54j . the corresponding quantum phase transition was suggested to be of the 2nd-order. 



Recently, similar phase transition within the fermion model for identical nucleon system has 



also been performed 



45 



46| 



From the periodic chart, one can deduce that nuclei that display an S0(6) spectrum lie 
close to the end of the shell, at least in the neutron sector. Therefore, to explore whether the 
transitional patterns between vibration and 7-soft spectrum can be realized in the SDPSM, 
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tion. 



T 1 , 1 , 1 , 1 r 




FIG. 3: B(E2) ratios vs a in the vibration-rotation transition. 
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TABLE II: Energy and B(E2) ratios at vibrational, rotational limit, and X(5)-like critical point 
calculated in the SDPSM. 



limit 




^1 
^2+ 


"a 


4+^9+ 
2+^0+ 


2+^0+ 


vibrational limit 


1.99 


2.97 


1.47 


1.49 


1.48 


X(5)-like point 


2.91 


5.60 


1.05 


1.38 


1.38 


rotational limit 


3.33 


6.96 


0.46 


1.34 


1.32 




"2 






2+^0+ 


4+-.2+ 




^2+ 


^2+ 


^2+ 


2+->0+ 


2+->0+ 


X(5)-like point (0^ band) 


5.32 


2.30 


5.33 


0.37 


0.43 



we considered a system with A'^ = A'^y = 3 in the gds shelL Namely, neutron pairs in this 



49|. By fitting 



case were treated as three neutron-hole pairs and a negative k was used as in 
i?42 = 2 and 2.5 for vibrational and 7-soft limiting cases, the parameters were fixed, and 
the results are listed in Table [H The detailed discussion about the two limiting cases in the 
SDPSM can be found in 3- 

The IBM calculation show that the level crossing-repulsion behavior of 0^ and Oj]' 
occurs 551] in the critical region of the f/(5)-5'0(6) transition. The SDPSM results of 0^ 
and Ojj' states, given in FiglU show the similar behavior of level crossing-repulsion when 
a = 0.58. Therefore, to see the behavior of effective order parameters against the control 
parameters clearly, the quantities related to Oj" state were also calculated for the 0;^ state. 

The results for -R42(-R6o)5 Ki{K2) and V2{v2) are given in FigJSl Fig|6]and FigjTJ respec- 
tively. The effective charges were fixed as e^^ = —^e^ = 1.5e since the neutron pairs were 
treated as holes. 

Figin^ shows that the typical ratios, -R42 = 2 (when a = 0) and 2.47 (when a = 1), of 
vibration and 7-soft spectra were produced. Interestingly, we found that in comparison with 
that of the rotation-vibration transitional results, i?42 in the vibration-7-soft transitional 
region increases with a smoothly. 



From Fig|6]one can see that as the IBM results |40[| and i?42 given in Figj5^, the wiggling 
behavior in Ki is smoothed out in the vibration-7-soft transition. One can also see that 
because the structure of 0^ and O3 exchange at a ^ 0.58, the amplitudes of B{E2; 0^ — ^ 2f) 
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FIG. 4: Energy levels of 0^ and 0^ states vs a in the vibration-7-soft transitional region. 



and B(E2] 0'^ 2^) also exchange at this point. 



In 52|, the experimental data of Xe and Ba isotopes were analyzed, in which for smaller 
neutron numbers, ^^^'^^^Ba and ^^^Xe, the 0^ state was taken in the Rqq if its B{E2) decay 
was consistent with a = N — 2. It was also shown thatjssl B{E2; 0^ 2f)/B{E2; 2f — »• 
of) = 0.07 for ^^^Pt, while it is 0.81 for ^^^Pt. By considering these results, the 0^ state 
were taken in the Rqo, K2 and V2 when a > 0.58. In comparison with those in the vibration- 
rotation transition, FigJSb and Figj6^ show that R^o, Ki, and K2 change smoothly with a, 



40 



52|. 



which are the typical features of the 2nd-order phase transition 

FigJT] shows that as predicted in the IBM and shell model calculation for identical system, 
the vibration-7-soft phase transition takes place and it is the second order phase transition, 
for which V2 and V2 change smoothly with a, the wiggling behavior changing sign in the 
region of the critical point are smoothed out. 



In the U{5)-SO{6) transitional region in the IBM, E{5) is the critical point symmetry 54, 
561 ]. It is interesting to see whether the signature of the -E(5)-like symmetry can be realized 
in the SDPSM for proton-neutron coupled system. We found that E^+/E2+ = 2.19 when 
a = 0.54 corresponding to the typical value of E{5) symmetry in the IBM. Other typical 
results are listed in Table UTTl in which the IBM results for = 5 are also given 5J]. It 
is seen that except for EQ+/E2+ = 2.59, which is smaller than that of the IBM result for 
= 5, the properties of -^(5) symmetry in the IBM indeed occurs in the SDPSM. 
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TABLE III: The SDPSM results for £^(5)-like symmetry. The corresponding results with = 5 in 
the IBM are also given 54 ] . 



limit 


^4+/-^2+ 


-^o+/-^o+ 




SDPSM 


2.19 


0.99 


2.59 


IBM 


2.19 


1.04 


3.68 




2+^0+ 


2+^2+ 
2+^0+ 


0+^2+ 
2+^0+ 


SDPSM 


1.36 


1.29 


0.53 


IBM 


1.38 


1.39 


0.51 




0+^2+ 
0+-2+ 


0+^2+ 
0+-2+ 




SDPSM 


0.06 


0.03 




IBM 











IV. SUMMARY 

In summary, the shape phase transition patterns for proton-neutron coupled system were 
studied within the framework of the SD-pair shell model. The results show that patterns 
of vibration-rotation and vibration-7-soft shape phase transitions are indeed similar to the 
corresponding results obtained from the IBM previously. The signatures of the critical point 
symmetry in the SD-pair shell model are also close to those shown in the IBM. The procedure 
may be extended to study quantum phase transitions in other fermion systems. 

This work was supported in part by the Natural Science Foundation of China (10675063; 
10775064), the U.S. National Science Foundation (0140300; 0500291), the Education De- 
partment of Liaoning Province (20060464), and the LSU-LNNU joint research program 
(LSU-9961). 
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